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Abstract 

In this paper we investigate the evolution of the viscous dark energy (DE) interacting with the 
\ dark matter (DM) in the Einstein cosmology model. Using the linearizing theory of the dynamical 
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Q ' system, we find, in our model, there exists a stable late time scaling solution which corresponds to 

the accelerating universe, and we also find the unstable solution under some appropriate parame- 
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ters. In order to alleviate the coincidence problem, some authors considered the effect of quantum 



account the bulk viscosity of the cosmic fluid, the coincidence problem will be softened just like 



\ correction due to the conform anomaly and the interacting dark energy model. But if we take into 
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. the interacting dark energy cosmology model. That's to say, both the non-perfect fluid model and 

\ the interacting models of the dark energy can alleviate or soften the singularity of the universe. 
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I. INTRODUCTION 

It's well known that recent data from la supernova (SN la) [l| and microwave background 
(CMB) radiation have provided strong evidences for a spatially fiat and accelerating uni- 
verse in the present time. The origin of accelerating expansion is regarded that the universe 
is dominated by an exotic component with the negative pressure called "dark energy" which 
constitutes 70 percent of the energy density of the universe, and dark matter about 26 per- 
cent. There are several candidates for dark energy: The first is the cosmological constant 
3], and the second is the so-called dynamic candidates such as: Phantom quintessence 
5], K-essence {g] and quintom The difference of these candidates for dark energy is the 
size of the parameter ue, namely the ration of the pressure and energy density of the dark 
energy. For quintessence, the state equation is given by the relation between the pressure 
Pe and the energy density pe, i.e. pe = ^^ePe, where — 1 < u;^; < —1/3. The borderline 
case of a;^; = — 1 of the extraordinary quintessence covers the cosmological constant term. 

The negative pressure of the dark energy may be the cause of the acceleration of the 
present Universe. However, the nature of the dark energy still remains a complete mystery. 
No more than five years ago, some physicists found that, if we only assumed the cosmic 
fluid to be ideal, i.e. nonviscous, it must bring out the occurrence of a singularity of the 
universe in the far future. There are two methods to modify or soften the singularity. The 
first is the effect of quantum corrections due to the conformal anomaly The other is to 
consider the bulk viscosity of the cosmic fluid flo|. The viscosity theory of relativistic fluids 
was first suggested by Eckart, Landau and Lifshitz In recent years some physicists 



12l | also took into account the bulk viscous cosmology. In this paper we consider that 



the bulk viscous dark energy is characterized by energy density pe and pressure as 
Pe = ujePe + PPe- 

Up to now the model, which the dark energy is considered as the perfect fluid state, 
suffers the coincidence problem [l^. That's to say, why the dark energy (DE) and the dark 
matter (DM) are comparable in size exactly right now. In recent years some interacting 
DE models have been brought out to overcome the problem. In these models they all are 
assumed that there exists a nonzero interaction between DE and DM in the Universe and 
gauges DE transfers to DM which allows us to create an equilibrium balance in the evolution 
of the Universe, so that the density of DE keeps the same order as that of DM at late times. 
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Some authors [ij] investigated dynamical behaviors of the dark energy models with only 
the dark energy linear equation of state interacting with dark matter in different cosmology 
models. They all found that the universe will enter an era dominated by dark energy and 
dark matter with interaction between them, and accelerate in late time under some proper 
parameters of the dynamical system. In this paper we focus on extending the equation of 
state of dark energy to nonlinear term, i.e. bulk viscous fluid, to investigate in what way 
the nonlinear term affects the evolution of the cosmology. 



II. DYNAMICS OF THE INTERACTING VISCOUS DARK ENERGY MODEL 
IN THE FLAT FRW UNIVERSE 

Base on the Einstein General Relativity theory, the standard Friedman equation acts as 
follows in the fiat FRW universe 

H'=^-^{Pm + Pe), (1) 

and the conservation relation of the cosmological total energy is 

p + 3H{p + p) = 0, (2) 

where G is the gravitational constant, and the the total energy density is p = Pm + Pe- 

From a hydromechanical standpoint a generalization of cosmic theory so as to encompass 
viscosity is most nature. In recent years Misner et.al [lOi] interestingly investigated the 
viscous cosmology. In this paper we consider the general equation of state equation of the 
viscous dark energy as 

Pe = ujePe + 13 p%- (3) 

Under considering the interaction term Q between the dark energy and dark matter, the 
evolution equations are 

Pe + 3if[(l + uje)pe + I3p%] = -Q, (4) 

pM + 3HpM = Q- (5) 

Differentiating Eq.(l), then putting Eq.(2) into it, we can get 

H = — {pM + Pe+Pe)- (6) 



If we introduce the follow dimensionless variables 

SnGpE SnGpM did . , 

where A'^ = Ina is the number of e-folding to present the cosmological time. The interacting 
term Q = 3bHpE is between the dark energy and dark matter. When the the coupling 
constant b is positive, it means that the dark energy converts into dark matter. Under these 
conditions we can construct the follow autonomous dynamics system by Eqs.(l), (3), (4) 
and (5) for the interacting viscous dark energy model in the flat FRW universe. 



x'^f{x,y) = -?>{l + b + ujE)x- (5x'^ 

+x[3y + 3(1 + uje)x + (8) 
y' = 9{X: y) = Sbx -3y + y[3y 

+3(1 +u;e)x + (3x% (9) 

The dynamics is the general form X' = F{X), where X is the column matrix constituted 
by the auxiliary variables and the prime denotes derivative with respect to = Ina. In 
order to analyze the stability of the dynamics system, we must linearly expand the dynamical 
system near the critical points due to the linearizing theory of dynamical system. So we can 
acquire the stability properties of the dynamical system from the eigenvalues of linearizing 
matrix. 

In the following section wc will focus on evolution of the dynamical system in the phase 
space and analyze its stability. 

III. EVOLUTION AND STABILITY ANALYSIS OF THE AUTONOMOUS DY- 
NAMICAL SYSTEM 

In order to investigate the evolution and its stability of the given dynamical system in the 
phase space, we take the parameter d = 2 with lowest order nonlinear term of the viscous 
dark energy Pe to its density pe for simplicity. We firstly solve the equations x' = f{x,y) = 
and y' — g{x,y) — to get the critical points: 

Point I: (x^,y^) = (0,l), 

and Point II: (x^J VJ) = ^^^), 
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where 7 = ^ ii"^ + 12/?fe + ^^u)e + 9a;|. 

For the Point I, it means that the universe is dominated by the dark matter. The 
eigenvalues of the hnearizing matrix for Point I: 

A( = 3,A^ = -3(6 + cjs). (10) 

We can see that when h < —uje, the two eigenvalues are all positive, so the critical point 
I is saddle. When b > —uje, so the critical point I is node. So we find Point I is an unstable 
critical, that's to say, the universe only dominated by the dark matter is unstable. 

For the Point II, the existence of the critical point II denotes to the era dominated both the 
dark energy and dark matter in the late-times of our universe, which is the same as classical 



Einstein cosmology [15|]. In Fig.l we showed the phase diagram of the autonomous dynamical 
system in the (x, y) phase space with the parameters P = 0.5, b = 0.8, ue = —1.2, d = 2. 

In order to investigate the stability of the critical point II, we linearize the dynamical 
system near the critical point II as follows: 

Sx' = [-3{b+l + ujE) + ^yc + H^ + ujE)xc 

-pdxp^ + (3{d + l)x^](5x + 3xc5y, (11) 
5y' = [3b + 2{l+UE)yc + (3dycx''c-^]5x 

+ [-?> + Qyc + 'i{l + u:E)xc + fix''c\5y. (12) 

From the above linearized equations, we can define the following linearizing matrix 

/ ^ I ^1 \ 
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dx 


\xc,yc 


dy 


\xc,yc 1 


1 dg 




dg 




\ dx 


\xc,yc 


dy 


\xc,yc j 



where the four elements of the matrix are 



9/ I 

dx '"^'^^ 
^1 

dy '"^'^^ 
dg , 
'"'''^^ 

dg 

d^ '"^'^^ 



-3(6 +1 + uje) + 3yc + 6(1 + uje)xc 

-(3dxf7^ + (3{d+l)x'^c^ 
3xc, 

36 + 3(1 + u;E)yc + Pdycx'}7\ 
-3 + 6yc + 3(1 + uje)xc + (3x'^c- 



(14) 
(15) 

(16) 
(17) 




FIG. 1: The phase diagram of the interacting viscous dark energy model in the flat FRW universe 
with the parameters (5 = 0.5,6 = 0.8, loe = — 1.2, d = 2. The point Cu is the critical point of the 
dynamical system. 
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FIG. 2: Critical point coordinates (a;c(solid) , yc'(dash)) for fixed /3 = 0.5,6 = 0.8 and d = 2. 
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FIG. 3: Critical point coordinates (.(('(solid) . yc'(dash)) for fixed 6 = 0.8, loe = —1-2 and d = 2. 
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FIG. 4: Critical point coordinates (xc(solid) , yc'(dash)) for fixed /3 = 0.5, oje = 



-1.2 and d = 2. 



The eigenvalues for the hnearizing matrix of the dynamical system near the critical point 
II are: 

Af = 2^(6/? + 2p^ + 18pb + 9pujE + 9cj| - 2^7 + 

-(36/32 _ + 3Qi3%^ _ 3Q/3^uje + 3Q(3%uoe 

+ IS/^^cul + 2lQ(3hujl + 108/?u;| + lQ2uj% 
-m^^ujE + 3QplriujE + 18/37cj| + 547u;|)^), (18) 

A^^ = ^(6/? + 2/?2 + i8/?6 + 9/3cJs + 9cj| - 2/?7 + ?>-iU}E 
+ (36/?2 - 72/?26 + 36/3%'^ - 3Qf3'^ujE + 3Qf3%ujE 
-lOSpul + ISP^ul + 216/36a;| + lOSPtol + 162u% 
-36p-fuJE + 36pb-fuJE + ISp-fLul + 547a;|)5). (19) 




FIG. 5: Eigenvalues A( (left column) and (right column) of linearizing matrix of the dynamical 
system near critical point II for fixed uje = —1.2 and d = 2. The second row is viewed through the 
h axial direction on the first row picture. 

Because of the complexity of the above two eigenvalue expresses, we can not simply 
determine them positive or negative depending on the parameters. In this paper we mainly 
perform the numerical simulations of Eqs.(18) and (19) in Fig. 5. From Fig. 5 we can see that 
there are same region of parameters /?, h where the eigenvalues Ai and A2 are both negative. 
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This means that critical point II is stable point as showed in Fig.l. Moreover if we investigate 
the stable point position (xc, yc) with the parameters which showed in Figs. 2-4, we can find 
that the universe will be dominated by dark matter both the interaction coefficient b and 
the viscous fluid coeflicient (3 become stronger. At the same time, for the critical point II, 
we can see that the coordinates {xc,yc) in the phase space is not vanished in its stable 
region, which tell us that the coincidence problem will be alleviated in the universe, and we 
also can see that the interaction term and the viscous dark energy will do the same effect 
on the alleviating the cosmological coincidence problem. 

IV. CONCLUSIONS 

In this paper we have investigated the evolution of the viscous cosmology model which 
dark energy interacts with dark matter. Using the linearizing theory of dynamical system, 
we found, in our model, there exists a stable late time scaling solution which corresponds to 
the accelerating universe, and we also found the unstable solution under some appropriate 
parameters. 

We all know that, in order to alleviate the coincidence problem, some authors considered 
the effect of quantum correction due to the conform anomaly: such as dynamical Casimir 
effect with conformal anomaly, or dark ffuid with conformal anomaly [9]. And some authors 
accounted some interacting dark energy models which can also soften the coincidence prob- 
lem [l^. In this paper that we found that, if we take into account the bulk viscosity of 
the cosmic ffuid, the viscosity softens the coincidence problem as the interacting dark en- 
ergy cosmology models. That's to say, both the non-perfect ffuid model and the interacting 
models of the dark energy can alleviate or soften the singularity of the universe. 
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